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Abstract 

We show that the large corrections due to final state interactions (FSI) in the 
TT~TT~^TT~^ , TT~TT~^TT~^ , and K~TT^'K^ decays can be accounted for 

by invoking scattering amplitudes in agreement with those derived from phase shifts 
studies. In this way, broad/overlapping resonances in S-waves are properly treated 
and the phase motions of the transition amplitudes are driven by the scattering 
matrix elements determined in many other experiments. This is an important step 
forward in resolving the puzzle of the FSI in these decays. We also discuss why 
the a and k resonances, hardly visible in scattering experiments, are much more 
prominent and clearly visible in these decays without destroying the agreement 
with the experimental tttt and Kir low energy S-wave phase shifts. 



T^^email: oller@uni.es 



1 Introduction 



During the last few years a series of theoretical works that combine the Chiral Perturbation 
Theory (CHPT) expansion [H |2j with unitarity in a consistent way EJ |31 13 IHl E] , 
order by order, have showed up the close relationship between the leading order CHPT 
amplitudes (driven by the chiral structure and the actual value of the order parameter for 
the spontaneous breakdown of chiral symmetry, /^r) and the presence, mass and width, 
together with other properties, of the lightest scalar resonances. These results, concerning 
the (J resonance, have been also confirmed by the solution of the Roy equations PU] • 

Also recently, but from an experimental point of view, charmed decays have revealed 
as a powerful source for the study of the properties of the scalar resonances offering 
experiments with high statistics, which is of foremost importance taking into account the 
traditionally lack of high statistics experiments in the scalar sector. In this respect, the 
studies of the Fermilab E791 Collaboration concerning the decays tt^tt+tt"'" [TT] . 

1686 events, and K^-k^tt^ ^21; 15090 events, were the first to point out statistically 

significant contributions of the a and k resonances, respectively. Other Collaborations 
on D decays also report with high statistical significance the existence of a cr resonance 
from the K'^tt+tt- decays, like the CLEO Col. p^, 5299 ± 73 events. Belle Col. 

[Hj, 57800 events, and the BaBar Col. [T3], 81396 events. Thus, we can conclude that 
nowadays, with the new data from charmed decays, the existence of the a should be taken 
for grant. 

Despite these positive facts concerning the charmed decays, from the theoretical point 
of view several questions arise regarding the analyses followed in the previous references, 
particularly in the S-waves. We concentrate in this work in the results of the E791 
Collaboration with respect to the 7r^7r+7r+ ^^1; 7r^7r+7r+ JT] and 

K^Ti^-n^ [12] decays. Although the final state is a three body one, we would expect 
that, at least, at the low energy tail of the invariant mass of the neutral 7r"'"7r~ subsystem, 
where the a is observed, the interactions of this two pion subsystem with the other pion 
should be quite soft due to the largeness of the mass. In this way, the movement of 
the phase of the a resonance across the Dalitz plot should be given, according to Watson 
final state theorem, by the isospin (/) 1/2 S-wave vrvr phase shifts. However, the Breit- 
Wigner (BW) employed for the a meson in ref.|,llj, as it is well known, does not fulfill this 
property, as we discuss below in more detail. Analogous comments are also appropriate 
for the K resonance in the K'tt^tt^ decay ^2] for the K'tt^ subsystem and the Kir 
7=1/2 S-wave phase shifts. In addition, controversial properties were also reported by 
this collaboration regarding the /o(980) and the i^o(1430) resonances. We will develop 
in this work alternative parameterizations, based on CHPT plus unitarity (chiral unitary 
approach [T7| UHl CH]), to those employed by the E791 Collaboration that do not show 
the mentioned problems with the S-waves but also contain the a and k resonance poles. 
We will also discuss, based on the presence of Adler zeros, why large destructive non- 
resonant contributions (called backgrounds by theorists) that completely distort the peak 
shapes of the a and k resonances are present in vrvr and Kti S-wave scattering at low 
energies, respectively. We will also offer a reason why these backgrounds are not present 
in charmed decays (also applicable to B decays) due to the lack of such zeroes occurring 
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only at a specific energy. As a result, we will put on sounder theoretical grounds the 
important findings of the E791 Collaboration about the relevant role played by the a 
and K mesons.*^ Let us stress that in this work we focus on the solution of the previous 
problems associated with the S-waves but we do not attempt to account for full three 
body unitarity, this is beyond the scope of this work. 

The content of this paper is as follows. After this brief introduction, we consider the 
— s> 7r^7r+7r+, 7r^7r+7r+ and K^ti^ti^ decays in sections |2l El and |^ 

respectively. We review the formalism and main findings of refs. |1H ITBl IT^ regarding 
these decays, in order, and develop our formalism to take into account FSI in section |21 
and apply it to the previous three decays. We also discuss why one should expect clear, 
although broad, peak structures associated with the a and k, resonances in these decays. 
We end with section El pointing out the most relevant results of our study. 

2 FSI in the D+ 7t-7t+7t+ Decay 

The study of the Dalitz plot of the decay — > 7r^7r+7r+, with a sample of 1686 candi- 
dates, was performed in ref.[llj within the Fermilab E791 Collaboration. The estimated 
signal to background ratio is 2:1, hence we are left with 1124 events and in the study of this 
decay we will normalize our results to this number of events. Previous studies with lower 
statistics can be also found in this reference from the E687 Collaboration. In order to in- 
troduce the controversial situation regarding the amplitude employed for describing data 
in ref. pi], let us briefiy discussed the fitting process followed by ref.^llj. In this reference, 
when the decay amplitude is modeled as the coherent sum of well established resonances 
pUj . following the isobar model, the non-resonant term is dominant while the description 
of the fit is poor. Indeed the psr degree of freedom, x^/^^ turns out to be 254 for 
162 degrees of freedom |ri|. The main discrepancies between data and the parametrized 
amplitude come, particularly, from the low energy region, below 0.5 GeV^. In order to 
improve the quality of the fit of the Dalitz plot the authors of ref. ^| included another 
tt^tt" resonance corresponding to the a resonance whose mass and width were allowed to 
fioat in the fit. The quality of the fit substantially improves with a x^/^ = 138/162 and 
finds values of 478 ± 24 and 324 ± 42 MeV for the mass and width of the a, respectively. 
However, the phase of the Breit-Wigner (BW) employed to describe the a resonance does 
not follow the elastic S-wave / = vrvr phase shifts. This is shown in Fig^ by the dashed 
line in comparison with the data points corresponding to the elastic S-wave / = tttt 
phase shifts from several experimental references. Indeed, the errors show the variation 
in the values from one reference to another, see ref. [Ij for more details. The discrepancy 
is manifest, specially close to threshold where the variation of the phase for the BW is 
much faster than that corresponding to data, despite that at the shown energy the vrvr 
interaction is elastic since the KK is above 1 GeV and the 4tt only starts contributing in 
a significant way above around 1.3-1.4 GeV. 

Denoting by 1 the n~ and by 2 and 3 the equal positively charged pions, the amplitude 

#2 Also called /o(600) and if* (800) in the PDG respectively. 
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employed in ref.^T] can be written as: 



N 



A = aoe'^'ATo + J] a„e*^" A(si2, si3)Ar„ . (2.1) 



n=l 



We describe in detail the ingredients present in the previous equation since we found 
several errata and omissions in ref.pi] that are necessary in order to reproduce their 
amplitudes. In Eq. ()2.1|) the first term is the non- resonant contribution and the other 
ones originate from the exchange of resonances. Every resonant contribution is Bose 
symmetrized for the equally charged pions, An = ^n[(12)3] +v4n[(13)2], as usual. The 
parenthesis around 12 mean that the particles 1 and 2 form the resonant state, and 
analogously for (13). The coefficients and 5„, for n > 0, are real constants that float in 
the fit, the a„ are magnitudes and the 6n phases. Finally, the A/^, n > 0, are normalization 
factors*^ given by: 

Xn = 1/ Q dSi2dSi3\An{Sl2,Si3)\^^ . (2.2) 

The ^„[(12)3](si2, S13) amplitudes for n > 1 correspond to the product 

An[{l2)3]{su,sis) = BWn{su)Fi,'\su)Fy\su)Ml'\{12)3] . (2.3) 
The Breit-Wigner propagator is given by, 

BWnisu) = {si2 - ml + imnTn{s 12)) ^ , 

r„(si2) = Tn{mJ^= ( — ] ^(j) , (2.4) 



S12 \PiJ F^'iPi) 

with rrin and r„(si2) the mass and 'energy dependent' width of the (12) resonance and 
Pi is the three-momentum of the particle 1 in the (12) rest frame. Furthermore, pi is pi 
evaluated at the resonance mass. On the other hand, the Aiif^ are angular dependent 
factors and read: 

M°)[(12)3] = 1, 
Alii)[(12)3] = -2pip3Cos^^i3 , 

M^^\12)3] = l(pip3)2(3cos^^i3-l) . (2.5) 

where pi and ps are the three-momenta of particles 1 and 3 and ^13 their relative angle, all 
of them referred to the (12) rest frame. The F^^ and are Blatt-Weisskopf penetration 
factors that depend on the spin J of the resonant state, and are given by: 



_ 1 
— I 



(0) 
D 

F^n^ = l/Va + Ns)^) 



'^'^ Omitted in the description of the amphtudes given in ref.|ll). but used in their results. 
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Fj^^ = l/V(9 + 3(gp3)2 + (rpi)4) , 

ir(o) = 1 

F« = l/v/(l + , 

= l/V(9 + 3(rpi)2 + (rpi)4) , (2.6) 

with g3 the three-momentum of particle 3 in the rest frame. 

The set of resonances that exchange in Eq. ()2.ip contain the p°(770), /o(980), /2(1270), 
/o(1370), p°(1450) and, in the parameterization that reproduces faithfully the experimen- 
tal data, the a. As discussed above the most controversial aspects present in Eq. (j2.H) 
involves the / = S-wave vrvr partial wave. 




0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.( 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

GeV GeV 
a) b) 



Figure 1: S-wave / = vrvr phase shifts (left panel) and modulus square of this partial wave 
(right panel), normalized such that the residue of the a pole is one. The data points are the 
elastic S-wave / = vrvr phase shifts from several experimental references and the errors show 
the variation in the values from one set to another, as explained in detail in ref.l^. The dashed 
lines correspond to the BW of the a resonance employed by the E791 Collaboration jllj . The 
thick solid lines are the full results of ref.^ when keeping only the vrvr channel while the thin 
ones correspond to the pure a pole contribution of Eq. (|2.7j) . The phase shifts and absolute value 
from the contribution of the a pole plus the first non-resonant term in Ea. H4.2() are shown by 
the dotted lines. Finally, the dashed-dotted lines are the results from the Laurent expansion of 
Ea. ()2.7|) keeping all the terms shown. 

Now, we want to show that we are able to reproduce the amplitude given in Eq. ()2.ip 
as employed in ref. but using the S-wave 1 = vrvr, KK coupled channel partial 
waves derived in ref.jnj- These T-matrices were obtained from Chiral Perturbation Theory 
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(CHPT) at leading order ^1 12] together with a unitarization scheme compatible with the 
chiral expansion (chiral unitariy approach). Furthermore, these matrices are not only able 
to reproduce the scattering data of the / = 0, 1 S-wave amplitudes up to around 1.2 GeV 
but also have been successfully tested by now in a vast number of production processes 
that pick up large corrections by FSI from these partial waves, see e.g. j^H 123 123 121j • 

The S-wave / = T-matrices from ref . [3| contains two poles in the appropriate Riemann 
sheets corresponding to the a and /o(980) resonances. Their pole positions are around 
448— i 224 and 990— ^13 MeV, respectively. Let us remark that the given poles are obtained 
in the full T-matrices derived in ref. j^, which are not any sum of pole contributions. These 
resonances originate because of the self interactions between the pseudoscalars and are of 
dynamical origin, disappearing in the large Nc limit jH |2H] • This is why they are said to 
be of dynamical origin. 

Due to the already mentioned discrepancy between phase shifts and the phase motion 
in energy of the BW for the o", let us perform a Laurent series around the a pole position 
of the / = TTTT S-wave amphtude, 

til = + 71 + l2{s -s„) + ... . (2.7) 

s- 

with the pole position and residua obtained from ref. [H] in the elastic case, with the values: 

= (0.47 - GeV^ , 7^ = 5.3 + i7.7 GeV^ , 

7i = -8.1 + i36.9 , 72 = 1.1 + i0.1 GeV-2 . (2.8) 

In Fig^we show in the left and right panels the phase and normalized absolute value of 
this partial wave, respectively. The thinner solid lines correspond to the pole contribution 
in Eq. ()2.7p . We see in FigH^ that the phase of this pole contribution does not vanish at 
threshold, but that runs parallel to the experimental phase shifts, so that the difference 
with respect to them keeps constant along energy up to the increase in the last points due 
to the closeness of the /o(980) and the opening of the KK channel. Thus, the phase of the 
cr pole contribution from Eq. (j2.7|) does follow the motion of the experimental S-wave 1 = 
TTTT phase shifts, in contrast with the BW phase indicated by the dashed line employed in 
ref.pi]. It is also worth stressing that the phase of the pure a pole contribution starts at 
—58 degrees at threshold and this is the reason why its value is not +90 degrees at the 
mass for the a, 466 MeV, but happens much later, around 1 GeV. By the same reason, 
for s 00 one gets only +122 degrees from this pole contribution instead of the usual 
+180 degrees. The agreement between the experimental phase shifts and the expansion 
()2.7|) is reached rather fast. The dotted lines in Fig|T] correspond to keep as well the 71 
term in Eq. ()2.7p . while the dashed-dotted ones correspond to keep all the terms shown 
in this equation. The thick solid lines are the full results from ref. [H] in the elastic case, 
removing the KK channel. In FigO^ we consider the absolute value of the partial wave 
under consideration, but normalized such that the residue at the pole position is one, 
so that the series in Eq. ()2.7|) is divided by 7g. The source of each line is the same as 
already explained and let us notice that this figure starts at GeV, below threshold. 
We want to stress three important facts: i) The markedly different behaviour of the a 
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Resonance 



Fraction 



NR 

cr7r+ 

pO(770)7r+ 

/o(980)7r+ 

/o(1370)7^^ 

/2(1270)7r^ 

p°(1450)7^^ 



(radians) 



0.57 
0.10 
1 (fixed) 
0.47 
0.24 
0.79 
0.20 



0.30 
3.40 
(fixed) 
2.90 
2.09 
1.05 
5.30 



3/152 



12% 
79% 
35% 
8% 
2% 
22% 
1% 



Table 1: Results of the reproduction of the parameterization of ref.^Jj summarized in Ea. ()2.1() . 
removing the BW of the a for its pole contribution, Eq. (|2.9jl . For every resonance we list the 
resulting magnitude a„ in the second column, the relative phase 5n in radians in the third column 
and the fraction of this decay mode in the fourth one. 



BW employed by the E791 Collaboration and the absolute value from the partial wave 
of ref.0. ii) The first peak in the BW. This peaks results because the BW formula ()2.4j] 
for the a resonance generates an unphysical pole in the physical sheet below threshold 
corresponding to a bound state at around 0.2 GeV. iii) The non-resonant contribution 
from the terms proportional to 71 and 72 in Eq. ()2.7|) are as big as the pole contribution and 
we see that the final shape of the absolute value of the amplitude is completely distorted 
as compared with the pure pole contribution. Indeed, the full results from ref. [S] show 
a zero at around 98 MeV. This is the Adler zero due to chiral symmetry, so that in the 
chiral limit the pseudoscalar interactions vanish at s = 0. At leading order in CHPT |2], 
this Adler zero sits as well at 98 MeV, very close to the position of the dip in the dotted 
and dashed-dotted lines and in good agreement with the zero in the total amplitude. In 
fact, because of the presence of this Adler zero, one can understand why the background 
turns out to be so big. If there is a pole that affects pretty much the low energy region, 
as in our case with the a resonance, it is then necessary a large background to cancel the 
pole contribution so that an Adler zero can happen. 

Before making use of the full results of ref.jS] let us substitute in Eq. ()2.H) the BW 
contribution given in Eg. (j2.4jl . as employed in ref. [TTj, by the pure a pole contribution, 

^ . (2.9) 

S - Sa 

located in the position given in Eq. ()2.8|) . This pole contribution, as discussed above and 
shown in Fig^, has a phase motion in agreement with that from the S-wave / = phase 
shifts. We keep the rest of terms in Eq. ()2.1|l and fit the Oj and 5i so as to reproduce the 
results from the parameterization employed in E791 with the values for the parameters 
given by their fit with the a. As in ref. ^T] the magnitude and phase of the p(770) vector 
resonance parameters, a„ and respectively, are fixed. In that fit apart from the a, and 
5j, the authors also leave as free parameters the mass and width of the a. Notice that we 
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do not include any Blatt-Weisskopf factors in Eq. ()2.9p . In order to reproduce the results 
from ref.jTTj we fit a Dalitz plot with 20 x 20 bins (a figure very standard in this E791 
analysis) normalized to the total number of events once the background is subtracted, 
namely 1124 events. This Dalitz plot is generated from the parameterization employed in 
ref.^T] for the signal only, corresponding to Eq. ()2.1|) . The resulting fit is very good with 
a low x^/^ = 3/152 and the obtained magnitudes and phases are given in Tabled On 
the other hand, we show in Fig|21the Si2 projection by the dashed line, while the results 
from the parameterization of ref. [TJj correspond to the points. We then conclude that we 
are able to reproduce the results of the E791 while keeping the constraint that the phase 
of the a contribution follows the / = S-wave vrvr phase shifts in a rather straightforward 
manner. An important difference of our fit in Table ^ with respect to the second fit of 
PI] is the clear dominance of the a pole in our case, 89%, as compared with the previous 
reference where its fraction, though also dominant, is lower, around 1/2. 



> 

CD 

o 
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Figure 2: m?{Tnr) projections with "data" points from Eq. (|2.H) of ref.'ll]. Our results are given 
by the dashed line with the pole contribution, Ea. H2.1j) . instead of the BW for the a of ref.|llj. 
The solid line corresponds to treat the S-wave FSI following Eq. (|2.13|) . 

Now, let us consider the full results of ref.p]. In this reference the partial waves are 
written as a product of two matrices (for the general case of coupled channels) as 

T = [I + N ■ g]~^ ■ N , (2.10) 

where tttt is channel 1 and KK is channel 2, in such a way that e.g. Tn is the 1 = elastic 
TTTT S-wave amplitude and so on. The diagonal matrix g{s) corresponds to the unitarity 
bubbles of every channel and the matrix N is determined by matching the previous 
expression with the chiral series at leading order, although the structure of Eq. ()2.10|) is 
valid to any order. At lowest order is the matrix of leading order CHPT amplitudes [3] . 
Now the important point for us is that given a general vertex Nij, projected in a certain 
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partial wave and connecting channels i and j, the summation of the unitarity bubbles just 
implies to multiply the matrix N by the inverse of the matrix [I + N ■ g], as follows from 
Eq. (j2.1(Jj) . In such a way that if the final channel i is produced from the initial channel k 
we then have: 

Tek = J2[I + N-9]jlN,k, (2.11) 

j 

so that j runs over all the possible intermediate states, that are also produced from channel 
k with the original amplitude Njk- Now we want to study a completely analogous situation 
where channel i is produced from a given production process starting with the meson. 
Then, if call by C,j the original amplitude for producing channel j the resummation of the 
unitarity bubbles must be given in exactly the same way as in Eq. (|2.11| ).^^ so that: 

^,-.J2[I + N-g]-'^, . (2.12) 

j 

In ref.0 the tttt and KK channels were considered for the study of the meson-meson 
/ = 0, 1 S-waves. The a and /o(980) appear as poles in the second Riemann sheet of the 
resulting partial waves. Thus, instead of considering Eq. (j2.H) . that includes explicitly the 
BW's associated with these resonances [11], we will instead apply Eq. (j2.12| ) to correct by 
final state interactions. In this way if we denote by D = [I + N ■ g], we will have: 

[Di,\si2) + D^,\sn)] a^.e''- + [D^iisu) + D^.^s^s)] aKRe'''''' , (2.13) 

in return of the contribution for the a and /o(980). Notice that we do not include in 
()2.13|) any form factor through Blatt-Weisskopf terms Fn'' and F^^ for the two scalar 
resonances. The rest of contributions, NR plus the ones from the /o(1370) and vector 
and tensor resonances in Tableware the same as in ref.jTTj, although their associated 
parameters and 5„ are fitted again. As in the case before, we fit a Dalitz plot with 
20 X 20 bins generated from ref.^J and normalized to the same number of total events. 
The resulting /v = 2/152 is of the same good quality as before, indicating an accurate 
reproduction of the amplitudes of ref.^I]. In Table|21we show the values of the parameters 
that we have fitted and in Fig|21 we give by the solid line the energy projections of any 
neutral rni subsystem. 

As a result we see that we are able to reproduce the signal function of ref. ^T] rather 
accurately and at the same time being able to establish that the scattering amplitudes 
which drive the final state interactions corrections in 7r^7r+7r+ are the same as those 

determined from scattering data and other production processes. The same conclusion 
is obtained in ref.^H]. This reference performs a K- matrix fit to data from and 
decays to Stt, although the resulting fits have significative lower confidence levels than 
those of the E791 Collaboration that we reproduce. Regarding the a resonance we see 
that it exhibits the same behaviour in the D-decays as in the rest of known processes, 
althogh in ref.j^H] is excluded. It is worth mentioning that the K-matrix employed in 

"^^The analytical properties of Nij and ^j, in general terms, are the same, that is, analytical functions 
except for the presence of cuts and poles. 
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Resonance 



Fraction 



(radians) 



NR 0.70 -0.43 17% 

(7r7r)7r+ 0.31 1.27 102% 

{KK)n+ 0.11 -0.39 6% 

p+(770)7r+ 1 (fixed) (fixed) 36% 

/o(1370)7r+ 0.31 1.99 3% 

/2(1270)7r+ 0.77 1.01 21% 

p°(1450)7r+ 0.20 5.48 1% 



2/152 



Table 2: Results of the reproduction of the parameterization of ref.^J summarized in Ea. (|2.1() . 
employing Eq. (|2.1H|l . For each intermediate state we list the resulting magnitude a„ in the 
second column, the relative phase 5n in radians in the third column and the fraction of this 
decay mode in the fourth one. 



ref.pi] does not meet the chiral requirement of a soft expansion for low energies and, in 
particular, it does not fulfill the chiral constraints imposed by the chiral power counting, 
which requires a vrvr scattering amplitude to start at order p^, instead of order p° as in 
This rises serious doubts about the applicatibility of the results from the K-matrix 
employed in ref . regarding the existence of the a resonance since meaningful stuctures 
are required in order to extrapolate the T-matrix in the energy complex plane away from 
the physical real axis, where it has been tested. 

Before ending this section let us discuss why the a resonance is more visible in the D 
decays than in vrvr scattering. One important point is the presence of the huge background 
(those terms proportional to 7i, 72, etc... in Eq. ()2.1|) ) in the / = S-wave vrvr partial wave 
as discussed at the beginning of this section. There we found out that this background is 
required in order to preserve the Adler zero in the presence of a light a resonance, so that 
a cancellation can occur close to threshold. In this way, the standard pole like structure of 
a resonance, e.g. that for the case the p resonance, is completely destroyed. The question 
now is why this background does not happen in D decays (or as well in B decays P^). 
The main point was already discussed in the second entry of ref. f2^, and arguments in 
this direction were also put forward in ref. |2Z] . The D ot B mesons can be identified with 
pseudoscalar sources coupled directly to a pion while the other two pions can be thought 
to couple just to a scalar source. As a direct application of CHPT power counting one 
realizes that this is not suppressed by any power of momentum or quark mass and then 
there is a priori no reason why the a meson should be screened by such large backgrounds 
as occur in the scattering case. Indeed, this is a result from Eq. ()2.13|1 . and can be seen by 
just making a Laurent expansion of the matrix around the a pole and then we check 
the absence of a significant background, contrarily to the scattering case. Certainly the 
scalar form factor of two pseudoscalars is not renormalization group invariant but this 
just amounts to a global quark mass multiplying factor and does not induce any energy 
dependence that could distort the pole structure of the resonance, contrarily to scattering 
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where the Adler zero happens for a specific energy. 



3 FSI in the Dj 77 77+77+ Decay 

In this section we consider the decay —>■ 7r~7r"'"7r"'" where the /o(980) resonance plays a 
central role as clearly seen in Fig|21 The "data" points correspond to the signal function 
of the E791 Collaboration given by the analogous formula to Eq. (j2.1|) but applied to this 
other decay. The set of intermediate resonances considered in ref.lj^ contains the /o(980), 
p°(770), /2(1270), /o(1370) and the p°(1450). In the fit performed in ref.UHl the mass 
and couplings of the /o(980) were allowed to float with an energy dependent width given 
by: 

^foi^) =9nVs/4:-ml + gk^ (^^ s/A~ m],^ + ^ s/A- m],}j , (3.1) 

to be substituted in the expression for BWnisu) in Eq. ()2.4|) . The striking fact is the 
small value for the coupling qk in the fit of ref.jTHI, compatible with zero, whereas the 
coupling of g.,^ is much larger. This result is very puzzled because is well known that the 
/o(980) has great affinity to couple with strangeness sources, as seen e.g. in decays 
|22j or in the J\l/ — > (pirn decay [23j. We show that we can reproduce the amplitude for 
the signal of ref. |16j making use of the amplitudes of ref. J3J, where the /o(980) has more 
standard properties regarding its couplings. We employ Eq. (j2.1|) of ref. [12] but, as in the 
previous section, we remove the /o(980) BW contribution of [12] and take into account 
the FSI by employing Eq. ()2.13p . in terms of the 1 = S-wave tttt and KK states. In 
order to reproduce the amplitudes of ref.[16j we proceed as in section |2] For that we fit 
a Dalitz plot with 20 x 20 bins, normalized to 625 events. The Dalitz plot is generated 
from Eq. ()2.H) with the parameters given in ref. JE]- We have performed two types of fits, 
given in Table (21 One in which the mass and width of the /o(1370) is fixed to the values 
of ref.lini and another in which, as in the previous reference, we let to fioat their values. 
The former corresponds to the second and third columns of Table (21 and the latter to the 
last two columns. Furthermore, in Fig|21we show the energy projection of any neutral tttt 
subsystem, the solid line refers to the former fit and the dashed line to the latter. 

We see that the reproduction of the signal function of ref . ^H] is very fair in both cases. 
In connection with the fore mentioned unexpected values of the couplings of the /o(980) 
as determined in ref.^lG^, we refer the interested reader to ref.[9j for a detailed study of 
the masses and couplings of the hghtest scalar resonances (a, k, /o(980) and ao(980)), 
showing that they obey rather accurately a standard SU(3) analysis and constitute the 
lightest scalar nonet. Other interesting studies on the /o(980) in and D decays are 
refs.[2EII2ni- 
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Resonance 


an 


8n 








Fraction 


(radians) 


Fraction(%) 


(radians) 


NR 


0.40 


0.16 


0.40 


-0.24 




13% 




14% 




(7r7r)7r+ 


0.28 


2.36 


0.25 


2.23 




6% 




5% 




{KK)tx+ 


1 (fixed) 


(fixed) 


1 (fixed) 


(fixed) 




78% 




84% 






0.24 


0.14 


0.22 


0.19 




4% 




4% 




/o(1370)7r+ 


0.60 


1.68 


0.57 


2.01 




28% 




28% 




/2(1270)7r+ 


0.50 


0.39 


0.48 


0.41 




20% 




19% 




p0(1450)7r+ 


0.25 


0.67 


0.24 


0.87 




5% 




5% 






11/142 


8.5/140 



Table 3: Results of the reproduction of the parameterization of ref. summarized in Eq. (|2.1|) 
for the Df tt^tt^tt^. For each intermediate state we list the resulting magnitude o„ in the 
second column, the relative phase Sn in radians in the third column and the fraction of this 
decay mode just below the value of every a„, analogously for the fourth and fifth columns. The 
second and third column correspond to the fit with the mass and width of the /o(1370) fixed to 
the values of ref. [16J . while in the fourth and fifth columns these values are allowed to float in 
the fit. In the latter case, nif^f^i^-jQ^ = 1.46 and ^ fg(i37o) = 0.16 GeV. 



4 FSI in the L>+ ^ K-n+n+ Decay 

The study of the high statistics Dahtz plot of the decay —>■ K^tt^tt^, with a sample 
of 15090 and an estimated background of a 6%, was performed in ref. [12j within the E791 
Collaboration in Fermilab. In order to introduce the controversial situation regarding the 
amplitude employed to describe data in ref.^l2j, let us briefly discussed the fitting process 
followed in this reference. When in this reference the amplitude, following the isobar 
model, is written as the coherent sum of well established resonances, as so quoted in the 
Particle Data Group ref. [20] , the description of the fit is poor. Indeed the per degree of 
freedom found in ref.^^l is x"^/^ = 167/63. The main discrepancies between data and the 
parametrized amplitude come, particularly, from the low energy region, below 0.6 GeV^, 
and at higher energies at around 2.5 GeV^. In order to improve the quality of the fit of 
the Dalitz plot the authors of ref. J2] allowed the mass and width of the scalar resonance 
i^'o(1430) to float in the fit and also included Gaussian scalar form factors a la Tornqvist, 
with the meson radii as additional free parameters. Despite that the x^/ is also around 
2. Finally, the authors of ref.[12j, together with the previous additions, also included 
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O 0.5 1 1 .5 „ 2 2.5 3 

GeV2 

Figure 3: mP^Tnr) projections for "data" ^B] and our results, solid and dashed lines, with the 
FSI given by Ea. H2.12() . In the dashed line the mass and width of the /o(1370) are allowed to 
float in the fit, second and third columns of Table 01 In the solid line the mass and width of the 
/o(1370) are fixed to the values of ref.flBj, fourth and fifth columns of Table |31 

the K resonance, referred as i^'o(800) in the last edition of the PDG 1201 tiut qualified 
as controversial. The quality of the fit substantially improves and the xV'^ = 46/63. 
However, the situation is far from being resolved since: i) the width of the well known 
Kq{14:30), allowed to float as discussed above, turned out to be a factor of two lower than 
its clearly determined value in scattering experiments and other production processes, ii) 
The Breit-Wigner employed to describe the k resonance does not follow the elastic S-wave 
1 = 1/2 Kir phase shifts. This is shown in FigE] where the points are the elastic S-wave 
1=1/2 Kir phase shifts from and the dashed line corresponds to the phase of the 
relativistic BW of the k. The discrepancy is manifest, specially close to threshold where 
the variation for the phase of the BW is much faster than that corresponding to data. 
Let us remember that the S-wave / = 1/2 Kit scattering is elastic below around 1.3 GeV, 
since both the Kt] and Ktitt channels are negligible below that energy, as clearly shown 
both experimentally [301 IHT| and theoretically [S], and then only the Ktt elastic channel 
matters in this energy region and the discrepancy cannot come due to inelastic effects. 

Denoting by 1 the K~ and by 2 and 3 the equal pions, the amplitudes employed in 
ref.p2] can be written as: 

N 

A = aoe''°Afo + ane''-An{su, s^^Wn ■ (4.1) 

n=l 

The different ingredients contained in the previous equation were discussed in detail in 
section |21 Here we only mention that in ref. the BW expression ()2.4|) has opposite 
sign. 
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The set of resonances that exchange in Eq. ()4.H) contains the K^{892), A'q(1430), 
ii'2(1430), i^'*(1680) and, in the parameterization that reproduces faithfully the experi- 
mental data, the k or i^'g(800). As discussed above the most controversial aspects present 
in Eq. ()4.1|) involves the J = 1/2 S-wave Kir partial wave. 

Now, we want to show that we are able to reproduce the amplitude given in Eq. ()4.ip 
as employed in ref. but using the S-wave 1 = 1/2 Kit, Kt] and Krj' coupled channel 
partial waves derived in ref. [S]. These T-matrices are obtained from Chiral Perturba- 
tion Theory (CHPT) |2I at next-to-leading order, supplemented with the exchange of 
explicit resonance fields in a chiral invariant manner [32], plus a unitarization scheme 
compatible with the chiral expansion (chiral unitariy approach). Furthermore, additional 
constrains from large QCD are considered as well in order to restrict the number of 
free parameters, for further details we refer to refs. jHl EHI- These T-matrices are able to 
provide an accurate reproduction of the Kti S-wave amplitude and of the J = 1/2 and 
7 = 3/2 S-wave phase shifts [101 EI] up to around 2 GeV. Later on, they were employed 
for calculating through dispersion relations the strangeness changing scalar form factors 
|33j . the light quark masses within QCD scalar sum rules [3^ and a crucial counterterm 
needed in the present and precise studies of K^^ decays [33] . 

The S-wave 1 = 1/2 T-matrices from ref. contains three poles in the appropriate 
Riemann sheets corresponding to the k, K^ilA'iQ) and the 7^0(1950). Their pole positions 
are around 708— i 305, 1450— i 142 and 1910 — z 27 MeV, respectively.*^ Let us concentrate 
on the first two resonances which are those that are relevant for the decay into 
K^Ti^-n-^, since the maximum value of the total center of mass energy of the subsystems 
(12) or (13) is mD+ — 171^,+ = 1.73 GeV, clearly below the infiuence of the ii'g(1950) pole. 
Let us remark that the given pole positions are obtained as a result of the full T-matrices 
derived in ref.|^, which are not just any sum of pole contributions, although these can 
dominate in some energy regions. The value for the width of the well established and 
clearly seen 7^0(1430) resonance is perfectly compatible with that quoted by the PDG 
[20], 294 ±23 MeV,#6 while that from the fit of the E791 Collaboration [12], 175 ±12 it 12 
MeV, is almost a factor of two lower. In addition, the mass obtained in ref.p2] for the 
7^0(1430) is out of the range given in the PDG as well. We also mention that the width 
of the K given by the E791 Collaboration, around 400 MeV, is substantially lower than 
the value from the pole position from ref . [H] , around 600 MeV. 

Due to the already mentioned discrepancy between the phase motion in energy of the 
BW for the k from the E791 Collaboration, let us perform a Laurent series around the k 
pole position of the 7=1/2 Kn S-wave amplitude, 

tn = + 71 + l2is -s,) + ... . (4.2) 
s- 

T^^These pole positions vary slightly depending of the fit taken from ref.jSj, although we have presented 
the values from the so called preferred fit, which behaves in the most appropriate way when considering 
scalar form factors and QCD sum rules. 

"i^^Remember that the double of minus the imaginary part of the pole position corresponds to the width 
of a resonance. This is indeed a possible and unambiguous definition of the width of a resonance. 
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Figure 4: S-wave 1 = 1/2 Kit phase shifts (left panel) and modulus square of this partial 
wave, normalized such that the residue at the k pole is one (right panel). The data points are 
from ref.pOj. The dashed lines correspond to the BW for the k resonance employed in ref.jl^. 
The solid lines correspond to the pure k pole contribution from Eq. (|4.2|) . The results from the 
contribution of the k, pole plus the first non-resonant term in Ea. H4.2|) are shown by the dotted 
lines. Finally, the dashed-dotted lines are the results from the Laurent expansion of Eq. (|4.2j) 
keeping all the terms shown. 

with the pole position and residua obtained from ref. |HI, with the values: 

= (0.71 - iO.Slf GeV^ , 7^ = 19.0 + i5.9 GeV^ , 
7i = 12.5 + i43.1 , 72 = 0.8 + i7.3 GeV"2 . (4.3) 

In FigE] we show in the left and right panels the phase and normalized absolute value 
of this partial wave, respectively. The thinner solid lines correspond to the pure pole 
contribution in Eg. ()4.2j) . We see in FigEK that the phase of this pole contribution does 
not vanish at threshold, but runs parallel to the experimental phase shifts, so that the 
difference with respect to them keeps constant along energy up to the increase in the 
last points due to the closeness of the Kq{14:30) and the opening of the Kr]' channel. 
Thus the phase of the pure pole contribution from Eq. ()4.2|) does follow the motion of the 
experimental S-wave 1 = 1/2 Kir phase shifts, in contrast with the BW phase indicated 
by the dashed line. It is also worth realizing that the phase of the pure k pole contribution 
starts at —90 degrees at threshold and this is the reason why its value is not +90 degrees 
at the mass of the k resonance, 708 MeV, but happens much later. By the same reason, 
for s —>■ 00 one gets only +90 degrees from this pole contribution. The agreement between 
the experimental phase shifts and the expansion ()4.2|1 is reached rather fast. The dotted 
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lines correspond to keep as well the 71 term while the dashed-dotted ones correspond to 
keep all the terms shown in this equation. The thick solid lines are the full results from 
ref.fS]. In FigEjD we consider the absolute value of the partial wave under consideration, 
but normahzed such that the residue at the pole position is one, so that we divide by 7q 
the series in Eq. ()4.2|) . The source of each line is the same as already explained and let 
us notice that this figure starts at around 0.3 GeV, below threshold. We want to stress 
three important facts: i) The very different behaviour of the k BW employed by the E791 
Collaboration and the absolute value from the partial wave of ref . jS] . ii) The first peak in 
the BW. This peaks results because the BW formula fl2.4|) for the k resonance generates 
an unphysical pole in the physical sheet below threshold at 0.46 ± iO.ll GeV, although it 
has a negligible impact above threshold, iii) The non-resonant contributions proportional 
to 7i and 72 in Eq. (j4.2j) are as big as the pole contribution, and we see that the final 
shape of the absolute value of the amplitude is completely distorted as compared with the 
pure pole contribution. Indeed they almost cancel each other at around 0.48 GeV giving 
rise to a zero in the full amplitude, as clearly seen in the figure. This is an Adler zero 
due to chiral symmetry, so that in the chiral limit the pseudoscalar interactions vanish at 
s = 0. At leading order in CHPT [2], this Adler zero sits at 0.48 GeV, very close to the 
position of the dip in the dotted and dashed-dotted lines and even closer to the zero in the 
total amplitude which also occurs at around 0.48 GeV. In fact, because of the presence 
of this Adler zero one can understand why the background turns out to be so big, in the 
same way as explained in section |21 for the a case. That is, if there is a pole that largely 
affects the low energy region, it is then necessary a large background to cancel the pole 
contribution so that the Adler zero can occur. 

Before making use of the full results of ref.0 let us substitute in Eq. (j4.H) the BW 
contribution given in Eq. ()2.4|) . as employed in ref. [12], by the k pole contribution. 



located in the position given in Eq. ()4.3|l . This pole contribution, as discussed above and 
shown in Fig|3^, has a phase motion in agreement with that from the S-wave 1 = 1/2 
phase shifts. We keep the rest of terms in Eq. ()4.1|) and fit the a, and 5i so as to reproduce 
the results from the parameterization employed in E791, Eq. ()4.1|) . with the values for the 
parameters given by their fit C. As in ref. the magnitude and phase of the i^'Q(892) 
vector resonance parameters, a„ and are fixed to 1 and 0, respectively. In that fit, 
apart from the Oj and the authors also leave as free parameters the masses and widths 
of the K and Kg (1430), as well as the meson radii that appear in the form factors. Notice 
that we do not include either any Blatt-Weisskopf factors nor Gaussian form factors in 
Eq. ()4.4|) . In order to reproduce the results from ref.^^l we fit a Dalitz plot with 20x20 
bins normalized to the total number of events once the background is subtracted, namely 
28369 events (the charge conjugate decay is also included). This Dalitz plot is generated 
from the parameterization employed in ref. JT] for the signal function, corresponding to 
Eq. ()4.H) . The resulting fit is very good with a low 'X^/v = 6.5/132. The values of the 
resulting magnitudes and phases are given in Tabled On the other hand, we show in Figl^l 
the S12 projection by the solid line, while the results from the parameterization of ref. [I^ 
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correspond to the points. The dashed hne is the so called high-projection (si2 > S13) and 
the dotted line corresponds to the low projection (si2 < S13). We then conclude that we 
are able to reproduce the results of the E791 at the same time that the phase of the k 
contribution follows the 7=1/2 S-wave Ktt phase shifts. 



Resonance 






Fraction 






(radians) 




NR 


2.10 


-5.95 


53.3% 




0.30 


-0.63 


56.7% 


K*(1430)7r+ 


1.00 


0.92 


12.2% 


i^i*(892)7r+ 


1 (fixed) 


(fixed) 


12.1% 


TT* (1430)71+ 


0.17 


-0.76 


0.4% 


ir* (1680)71+ 


0.45 


0.58 


2.5% 



X^/v 6.5/132 



Table 4: Results of the reproduction of the parameterization of ref.^2] summarized in Ea. ()4.1() . 
removing the BW of the k by its pole contribution, Eq. (|4.4j) . For each resonance we list the 
resulting magnitude in the second column, the relative phase 5n in radians in the third column 
and the fraction of this decay mode in the fourth one. 

We notice that in Tabled the so called NR contribution plays a more important role, 
it has a fraction of around 50% while in fit C of the E791 Collaboration amounts to just 
a 13%. 

Now, let us consider the full results of ref.|8] in order to take into account the FSI 
employing Eq. ()2.13p . In ref.0 the partial waves are written as a product of two matrices 
(for the general case of coupled channel) as T = [I + N ■ ■ N such that e.g. Tn is the 
1 = 1/2 elastic Krt S-wave. The diagonal g{s) matrix corresponds to the unitarity bubbles 
of each channel and the matrix is determined by matching the previous expression with 
the chiral series of CHPT plus resonances in the ?7(3) case, for more details on this respect 
see ref. [S^. In this reference the Kir, Kt] and Krj' channels are considered, hence we rewrite 
Eq. ()2.13p for this particular case as, 

+ [Diiis^^) + /^it(5i3)] ax,' e*'-"' . (4.5) 

Let us stress that, since the k and i^'o(1430) appear as poles in the second Riemann 
sheet of the partial waves of ref. |0, they are already taken into account when considering 
Eq. (j4.5|) . Notice that we do not include in Eq. (j4.5p for the two scalar resonances either 
any from factor as done in Eq. (j4.ip through the Blatt-Weisskopf terms and F^^ 
nor Gaussian form factors. Furthermore, the i^'o(1430) resonance from ref. |HI has a mass 
and a width in complete agreement with the values in the PDG [20], which are clearly 
determined from the study of Ktt scattering j30[ I31j. The rest of contributions present in 
Eq. ()4.ip , namely, the NR plus the ones from the vector and tensor resonances in Table |3] 
are kept, although the parameters and 5„ are fit again so as to reproduce the results 
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Figure 5: m?{K7r) projections for "data" ^ and our results with the pole contribution in 
Ea. ()4.1() instead of the BW for the k of ref.^^l- The solid, dashed and dotted lines correspond 
to our results for the projections nv^i^KTr^j jn^i^Kir^iQ^ and m'^{KiT)high, in order. 

of the signal parameterization of the E791 Collaboration. As in the cases before, we fit a 
Dalitz plot with 20 x 20 bins generated from ref . ^21 and normalized to the total number of 
events. The resulting x^/^ = 127/128 is acceptable, indicating a rather good reproduction 
of the amplitudes of ref.p2]- In Table El we show the values of the parameters that we 
have fitted and in FigElwe show the energy projections of any neutral Kit subsystem. 



Resonance a„ 5„ Fraction 

(radians) 

NR EeO (UO 29.6% 

{K7r)7i+ 1.66 4.10 31.8% 

lKri)7i+ 0.86 2.63 2.0% 

lKr]')TT+ 2.33 -1.54 9.8% 

Kl{892)7i+ 1 (fixed) (fixed) 11.6% 

ii'2*(1430)7r+ 0.11 -0.62 0.2% 

K*{1680)n+ 0.72 0.80 5.9% 

xViy 127/128 



Table 5: Results of the reproduction of the parameterization of ref. summarized in Ea. (|4.1() . 
employing Eq. ()4.5() . For each intermediate state we list the resulting magnitude in the second 
column, the relative phase 6n in radians in the third column and the fraction of this decay mode 
in the fourth one. For this fit we also leave as free parameters the poorly measured mass and 
width of the i^^(1680), their values are given in the text. 
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Figure 6: m?{KTT) projections for "data" ^ and our results taking into account Ea. (|4.5j) . The 
solid, dashed and dotted lines correspond to our results for the projections TvPi^Kir^, mPi^KTr^iQ^j 
and m?{KTr)high, in order. The dashed-dotted line refers to the fit when the {Kr])iT^ channel is 
excluded, Ea- lflT^ . 

As a result we see that we are able to reproduce the signal function of ref. J2] rather 
accurately and, at the same time, we are able to establish that the scattering amplitudes 
driving the final state interactions corrections in D'^ K~tt^tt^ are compatible with 
those determined from scattering data ^3^|M]. In particular, the k resonance is present 
in both and the mass and width of the i^'o(1430) does neither differ in both cases. On the 
other hand in the previous fit we have allowed to float the mass and width of the 7^^(1680), 
since they are poorly measured. In the PDG the values reported are between 1.7 to 
1.8 MeV for the mass and from 0.17 to 0.40 GeV for the width. The values that we obtain 
in the fit to the parameterization of ref . jT^ , which employs the central values given in the 
PDG, are m = 1.7 GeV and width T = 0.17 GeV. 

Although we have already included the Kt] channel in the fit of Table since this 
channel has little effect on the Kit scattering we then check the stability of this fit by 
removing the Kt] channel. Then, instead of Eq. (j2.13|) we now consider, 

[D^,\su) + (sis)] ax.e*'^- + [D^Hsu) + D^i{srs)] aK^,^'^^' . (4.6) 

The resulting x^/^ = 144/130 is a bit larger than the previous one, and the values of 
the parameters in the fit are rather similar so that we refrain from presenting them. 
Noneteless, we show the event projection in FiglHlfor this case by the dashed-dotted line. 
On the other hand, the resulting mass and width of the i^'*(1680) are 1.7 GeV and 0.13 
GeV, in order. Thus, we conclude that our results are rather stable under the presence or 
removal of the Krj channel although the resulting is somewhat lower when this channel 
is considered as well. 
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Finally, the same reasons advocated at the end of section |21 to explain why the a meson 
is clearly visible in decays to three pions can be also applied here for the k in the 

TT'TT+TT"'' decay.*'^ 

5 Conclusions 

We have considered in detail the FSI of the 7r^7r+7r+, —>■ 7r^7r+7r+ and 

— s> K^-K^n^ driven by the S-waves. The Dalitz plots associated with these decays 
were studied originally in refs. [TT| [W\ IT^ showing by the first time statistically signifi- 
cant evidences for the existence of the a and k resonances. Here we have payed special 
attention to those aspects still controversial after those works with the aim of improving 
the theoretical basis of the parameterizations employed in these references. In particular, 
we have shown that the meson-meson S-waves with / = 0, 1/2 that drive the correspond- 
ing FSI in the previous decays are compatible with those amplitudes determined from 
studies of scattering data [3 El |H], also tested in many other production processes, e.g. 
|211 122 I2ni 1211 ■ In particular, we have shown that the phase motion of the low energy 
FSI in the — >■ tt'tt^tt'^ and —>■ K~t:^tt^ decays follows the elastic / = tttt and 
1 = 1/2 Kti S-waves phase shifts, in order, and how this is compatible with the presence 
of the a and k resonances, respectively. We have seen that the reason for the disagree- 
ment between the phases associated with these resonances in the studies of refs. |1H IT^ 
and the experimental phase shifts is the employment by these authors of Breit-Wigner 
propagators. Once the BW propagator is substituted by the pure pole contributions from 
Eqs. (|2.7|) and (|4.2|) . for the a and k, respectively, the agreement is restored. Indeed, we 
have also shown that these pole contributions are not affected by significant backgrounds 
in the expansion of -D~^, see Eqs. ()2.13|) and ()4.5p . while huge destructive backgrounds 
are present in the scattering. In this way the pole contributions are not distorted in the 
studied Dg and D decays while they are so in the scattering, as shown in FigsH] and El 
Furthermore, we have also considered the full results for the S-wave FSI from refs. [21 El 
from Eqs. ()2.13|) and ()4.5|) . In this way we have considered simultaneously the interplay 
between the a and the /o(980) for / = (|2IIS1), and that between the k and ii'o(1430) 
for J = 1/2 ()4.5|) . This is also important since in ref. |TH] the properties of the couplings 
of the /o(980) are astonishing while in ref.^^ the mass of the i^o(1430) resonance is out 
of the range given in the PDG [201 and its width is almost a factor 1/2 smaller. Thus, 
our results show that one can also understand these FSI with "standard" /o(980) and 
7^0(1430) properties. 
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